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Abstract. - The statistics of signal increments are commonly used in order to test for possible 
intermittent properties in experimental or synthetic data. However, for signals with steep power 
spectra [i.e., E{uj) ~ a;~" with n > 3], the increments are poorly informative and the classical 
phenomenological relationship between the scaling exponents of the second-order structure func- 
tion and of the power spectrum does not hold. We show that in these conditions the relevant 
quantities to compute are the second or higher degree differences of the signal. Using this statis- 
tical framework to analyze a synthetic signal and experimental data of wave turbulence on a fluid 
surface, we accurately characterize intermittency of these data with steep power spectra. The 
general application of this methodology to study intermittency of experimental signals with steep 
power spectra is discussed. 
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Introduction. — Since the prediction of Kolmogorov 
in 1941 [T], it is well known that the spatial power spec- 
trum E{k) of a fluid particle velocity v in a turbulent 
flow is a power- law of the wave number k as k~^'^. 
The —5/3 exponent of the spatial power spectrum is re- 
lated to the second-order moment of velocity increments 
S2{r) = {[v{l + r) — w(/)]^) ~ r^/'^, I and r being a posi- 
tion and a spatial separation [5]. The phenomenological 
pelationship between both exponents comes from Fourier 
transform properties. It can be generalized to any station- 
ary random processes: if the power spectrum of the pro- 
cess is E{k) - fc-", then 52(r) ~ r^^ ^j^h (2 = n - I [2]. 
This property allows to perform measurements in the real 
space to reach the power-law exponent of the power spec- 
trum. The statistics of velocity increments are also cru- 
cial to characterize the intermittent nature of the velocity 
field using the scaling properties of structure functions: 
Sp{r) = {[v{l + r) — v{l)]P) ^ r^p {p positive integer) [3]. 
A non-linear dependence of (^p versus p is the hallmark of 
intermittency. 

Steep power-law spectra (^ aj~" or ^ fc^" with n close 
or larger than 3) of a process are observed in various 
situations: magnetohydrodynamics turbulence [?|, atmo- 



spherics turbulence [S], gravity [5] or capillary [7] wave 
turbulence on a fluid surface, and direct cascade of two- 
dimensional fluid turbulence [8l[9]. Whatever the corre- 
sponding signal measured in space or in time {e.g. fluid 
velocity or vorticity, surface- wave height, magnetic field, 
wind), such steep spectra mean that the measured sig- 
nal is at least once continuously differentiable [2]. The 
signal differences or increments are thus poorly informa- 
tive since they are dominated by the differentiable com- 
ponent of the signal. For instance, some numerical sim- 
ulations of the power-law scaling of the energy spectrum 
in two-dimensional turbulence exhibited apparent contra- 
dictions with its reconstruction from spatial correlation 
measurement, i.e. C2 7^ "- — 1 (see Ref. [10]). Babi- 
ano et al. have systematically reconsidered the theoretical 
relations between second-order structure functions and en- 
ergy spectra instead of phenomenological or dimensional 
arguments [TO]. They showed that the apparent contra- 
dictions come from the fact that the relation C2 = n — 1 
does not hold for steep power-law spectra. Indeed, due to 
the differentiable component, the exponent of the second- 
order structure function is independent of the spectrum 
slope as soon as this one is steeper than —3; that is C2 = 2 
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whatever n > 3 [TO]. This latter property has been noted 
elsewhere without derivation [5J[TTJ[T2]. The indirect con- 
clusions drawn from the structure function analysis to 
reach the exponent of the power-law spectrum (using the 
relation ^2 = n — 1) are thus misleading for ri > 3, but 
suprinsingly are still used for some experimental signals 
with steep power spectra |13) . 

In this letter, we emphasize that the increments are not 
the relevant quantities in order to statistically character- 
ize the fluctuations of a steep power spectrum signal and, 
in particular, to probe for its possible intermittent na- 
ture. We show analytically that, depending on the spec- 
tral steepness, it is necessary to adapt the degree of the 
difference statistics used to analyze the signal. Given an 
adapted degree, we provide a general relationship between 
the exponents of the second-order structure function and 
of the power spectrum whatever the steepness of the spec- 
trum. Finally, applying this approach to a synthetic sig- 
nal, and to an experimental signal of wave turbulence on a 
fluid surface allows us to accurately characterize intermit- 
tency of these data with steep power spectra. Note that a 
general framework for the study of intermittency of a sig- 
nal with arbitrary degree of regularity has been previously 
proposed using more complex estimators based on the con- 
tinuous wavelet transform [T3] or on inverse statistics [T^ . 
Here, we propose a practical approach and provide simple 
rules that should be easily applicable on the workbench 
to study intermittency of experimental signals with steep 
spectra. 

Scaling properties of irregular signals. Let us 

first recall the pioneering work of Parisi and Frisch [3] 
introduced for the description of the irregular nature of 
longitudinal velocity data in fully developed turbulence. 
They locally described the fluctuations of an erratic signal, 
77 (i), by means of the singularity exponents hito) which 
caracterizes the power-law behaviour of the finite differ- 
ences (or increments) of 77 over a time lag r at a time to 
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Since the characterization of non-local singularities cannot 
be achieved in a purely local manner, Parisi and Frisch 
introduced the p-order structure functions 



S,ir)^{\SMt)\'') , 



(2) 



where (•) represents an average over time t and from which 
they defined the spectra of global exponents C,p such as 



Sp{T) CX T<- 
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Note that it is useful to consider the moment of order p 
of |i5r'7(0| rather than 5rr]{t) as used in Ref. [3], so that 
Eqs. ([2]) and ([3]) are defined for all p positive real. In 
the presence of homogeneous fluctuations, i.e. h(to) = H 
whatever to, it is straightforward that ^p is a linear func- 
tion of p with (^p = pH. Reciprocally, a non-linear depen- 
dence of (^p with p is the signature of non-homogeneous 



fluctuations whose singularity exponents vary with time, 
i.e. the hallmark of intermittency. 

Scaling exponents of two point statistics. We 

now resume the relationships between the power-law expo- 
nents of the second-order structure functions, the correla- 
tion function, and the power spectrum. Let us consider a 
stochastics process 77(f) which is not necessarily stationary 
(see below). We assume that finite differences of 77 over 
a time lag r, (Sr'7(i) = vi^ + ''') ^ ^(^)i form a stationary 
process of zero mean and that 77(0) = 0. Let us look at the 
correlation between two short intervals (of size A) sepa- 
rated by a lag 9, that is the quantity (5A?/(i + ^)'^A?/(i)- We 
choose A equals to the sampling time of ri{t), that is the 
duration between two successive measurement points. For 
p ~ 2, the scaling of Eq. (JS)) implies a power-law decay of 
the correlation function C of the increments (5a 77 at large 
lags 9 as [lB| 

C{9) = {6AV{t + S)SAV{t)) « |6'r'^ ,withK = 2-C2 . 

|e|i-i.oo 

(4) 
This scaling behaviour for large lags coincides to a low- 
frequency power-law behaviour of the power spectrum C 
of the increments 5a 77, defined as the Fourier transform 
(FT) of the correlation function [TO] 



C{uj) = / C{e)e' 



M6' oc \uj\ 



with /3 = C2 - 1 , 



(5) 

u! being the angular frequency. An estimator of the power 
spectrum C{uj) is obtained taking the square modulus of 
the FT of the increments 5/^rj{t) observed over a finite 
time range [0,r], 



C(w) ~ Es^r^iuj) = 



5i,i']{t)e''^'At 



(6) 



In practice, a common habit is to compute the empirical 
power spectrum of 77, denoted E(lo), rather than that of 
the increments, E^^ni'^)- Indeed, these spectra are related 
by Es^,]{uj) = 2 [1 — cos(w)] E{uj) leading to the empirical 
power spectrum scaling 



E{lu) 
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2 [1 — cos(cl')] \^\^o 



oc |w|-" , withn = /3 + 2 . (7) 



The scaling exponents of S2{t) oc t^=^, C{9) oc \9\ " and 
C{uj) oc uj~" for the increments thus read respectively 



C2 



1, K = 3 — 71, and /3 



(8) 



where 77, < 3 (see below) is the exponent of the power 
spectrum of 'q{t) [see Eq. d?])]. 

It is fundamental to note that the above stationarity 
condition for the increments does not imply the station- 
arity of the signal, so that the correlation function of 77 
might not be defined and, thus, that E may not be a 
power spectrum in the statistical sense (i.e. the FT of 
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Table 1: Relationships between the spectral exponent n, the structure function exponent and the differentiability of a signal 



ri(t) with a self-similar power spectrum. Reference indicates the existing derivation of the relationships. ^2 (t) = {\5t 77] ) 
Sf (r) = {\S?^vf} and Sf (r) = {\Si'^v\'). 

Power spectrum DifFeren- Difference statistics used to test intermittency 
E{uj) - w~" tiability 



Second-order 
structure funct. 



n < 3 
n> 3 > 1 



71 < 5 1 
n> b > 2 



71 < 7 
71 > 7 



2 
>3 
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(1) 



S^^^T] = 7?(i + 2t) - 27;(t + t) + 7;(f) 
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a correlation function). Also, one should be careful that 
the estimation of the spectrum £'(0;) for a non-stationary 
signal r]{t) may be significantly biased depending on the 
FT numerical algorithm used. In such condition, a more 
robust practice to estimate the spectral scaling exponent 
n consists in numerically computing the FT of the sta- 
tionary signal (5A?7(i), i.e., Es^ri{,<^), and to use Eq. ([7]) to 
convert back to the usual power spectrum. 



Meaning of steep power spectra. Here, we ex- 
plain why the usual relationships of Eq. ([8|) no longer hold 
in the presence of steep power spectra (n > 3). In this 
case, Eq. (|S]) leads to k < 0, suggesting that the correla- 
tion function of SaV does not go to but rather diverge 
at large lag values [see Eq. Q]. Obviously, this is neither 
physically nor statistically acceptable. Indeed, for a non 
trivial stationary process (5,-f? satisfying Eq. ^, the corre- 
lation function of Eq. ^ is only defined for k = 2 — C2 > 0, 
that is for C2 < 2. This implies that the spectrum of the 
increments is only defined for /3 = C2 ^ 1 < 1- Con- 
sequently, for a process with stationary increments, the 
scaling exponent of the empirical power spectrum must 
satisfies n = C2 + 1 < 3. Hence, for steep power spectra 
(ri > 3), the basic assumption that increments form a sta- 
tionary process is not verified so that the structure func- 
tions of Eq. (j31) and the correlation function of Eq. (U) are 
not well defined. In practice, the classical phenomenologi- 



cal relation between E{u}) 



and ^2 (r) ~ r" ^ is thus 



invalid for n > 3, so that the spectral slope can not be de- 
duced from the measurement of the second-order structure 
function. This also means that the process r]{t) is at least 
once differentiable at times where //(i + r) — ri{t) ~ rdrj/dt 
at the first order in r [5]. These local linear trends are 
responsible for the non-stationarity of signal increments 
and, in turns, bias the estimation of scaling exponents. 
Indeed, near these times, one has \SAil{t)\^ ^ tP that cor- 
responds to (p = p [using Eqs. ^ and ^]. This means 
that the scaling of the exponent of the structure functions 
are independent of the spectral steepness. Thus, the in- 
crements of the signal do not appear as relevant quantities 
when looking for possible intermittency (i.e. a non-linear 
evolution of (^p with p) in signals with steep power spectra. 



Using higher-degree difference statistics to re- 
cover stationarity. — As recalled above, the scaling 
exponent of the spectrum of the increments is decreased 
by two with respect to the one of the spectrum of the sig- 
nal (/? = 71 — 2). Clearly, the repetition of the difference 
process allows recovering the power spectrum of a station- 
ary process. For instance, for 3 < 71 < 5, the second- 
degree difference of the signal S)^ ri{t) = SAlSAVit)] = 
r]{t + 2A) — 277(i + A) + r/(t) has a power spectrum with 
scaling exponent f3^'^' — n — 4 < I, which is compati- 
ble with S^^'T]{t) being stationary. The second-degree dif- 
ferences thus remove the local linear trends in the sig- 
nal 77 responsible for the saturation C,p = p when 71 > 3. 
Thus, when looking for possible intermittency in this case, 
one should use the structure functions of degree 2, that 



is S^p\t) 

A-2) 

'32 



{\5^r\{t)\P) OC 



['p . We have indeed 

(2) 



— 1 for 71 < 5. For ri > 5, one have Q 4 due to 
local quadratic trends in the signal. Thus, for 5 < ri < 7, 
the power spectrum of third-degree differences S^Jrjit) = 

5A[5^AV{t)] = V{t + 3A) - 37](i + 2A) + 3iit + A) - T]{t), 
is well defined and intermittency should be tested with 
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These results are sum- 



marized in the Table [T] Note that, even though the slope 
of the power spectrum is reduced by computing higher- 
degree difference statistics, the upper bound of the spec- 
tral bandwidth, related to the finite dynamic resolution of 
the original signal measurement, is not bypassed. 

Workbench recipe. — In practice, given a first esti- 
mate of the empirical power spectrum scaling exponent 71, 
further statistical analysis should be performed using dif- 
ference statistics of degree d > d* where d* is the smallest 
integer such that n — 2d* < 1. For instance, structure func- 
tion of degree d, Sp (r) = {\St 'r]{t)\P) oc |r|^p should 
be used to test for intermittency while spectral analysis 
or correlation analysis should be done on degree d differ- 
ences, Eg(d) (uj) oc \uj\^ and Cg(d) {9) ex \0\^'^ . The 
relationships between the scaling exponents then become 



Ad) 
'52 



1, 



K^''^ 



1 + 2d - n, /3('') = 71 - 2d, (9) 



for 71 < 1 + 2d. In other words, the proposed procedure 
removes the biases in the estimation of the scaling expo- 
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Frequency (Hz) 



200 



Fig. 1: Power spectra of the experimental data r]{t) (solid line) 
and of a synthetic signal (dash-dotted line). Dashed lines have 
slopes of -4.3 and -2.8. Inset: Typical temporal evolution of 
ri[t) during 10 s, (77) ~ 0. 



nents. These biases arise from the local regular behaviors 
of a signal that can occur at different degrees of the sig- 
nal differentiability. Hence, in order to numerically check 
that the signal regular components have been adequately 
removed, it is good practice to check that results remains 
consistent when increasing degree d to d+ 1. In particu- 
lar, one should make sure that Cp ~ Q . Note that 
for discontinuous signals the use of increments or higher- 
degree differences is unsuitable, a wavelet-based approach 
is more suited |14| . 

Applications. — To illustrate the results of the previ- 
ous section, let us now apply the proposed estimator based 
on higher-degree difference statistics to signals with steep 
spectra and probe their possible intermittent nature. A 
synthetic signal with prescribed intermittency and exper- 
imental data of wave turbulence on a fluid surface will 
be tested below for comparison. To our knowledge, only 
one study has compared the method of second-degree dif- 
ferences with more complex estimators based on inverse 
statistics in order to probe intermittency in a simulation 
of two-dimensional flows [17] . 

Synthetic data. Here, we apply the above suggested 
estimator to a synthetic data. Recently, it has been 
proposed that the scaling properties of experimental ve- 
locity (transverse to the mean flow) in fully developed 
turbulence could be described by log-normal Random 
Wavelet Cascade (RWC) [TS]. RWC generalizes the con- 
cept of self-similar cascades leading to multifractal mea- 
sures (— l<ri<l)to the construction of scale- invariant 
signals (n > 1) using orthonormal wavelet basis |18| . In- 
stead of redistributing the measure over sub- intervals with 
multiplicative weights, it allocates the wavelet coefficients 
in a multiplicative way on the dyadic grid. This method 
has been implemented to generate multifractal functions 
from a given deterministic or probalistic multiplicative 
process. From a mathematical point of view, the con- 
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Fig. 2: Rescaled second-order moment of the structure func- 



tions S2{t)/i 



with n = 4.3 computed from the (o) first- 



(D) second-, (o) third- and (A) fourth-degree differences of the 
signal as a function of the time lag r. a) Synthetic signal, b) 
Experimental signal. Correlation time is Tc — 63 ms. 

vergence of the cascade and the regularity properties of 
the so-obtained stochastic functions have been discussed 
in Ref . [T2] . Intermittency of RWC is characterized by the 
theoretical scaling exponents Cp = cip — C2P^/2 [TH]. Here, 
we consider a realization of 3 x 10® data points of the RWC 
process and choose ci = 1.92 and C2 = 0.27 to reproduce 
the intermittent properties of experimental data of wave 
turbulence (see below). 

We first compute the usual power spectrum of the RWC 
signals allowing us to assess that its frequency power law 
exponent is around 3 < n < 5. This suggests that unbi- 
ased estimates of the scaling exponents will be obtained 
using second (or higher) degree difference statistics (see 
Table [T]) . The unbiased power spectrum of this signal is 
then computed using the Fourier transform [see Eqs. © 
and (O] of the second-degree differences and is shown in 
Fig. [TJ It roughly behaves as a steep power-law over one 
decade in frequency, i.e. E{uj) (x u^^ with n ~ 4.3. The 
second-order structure functions S2 of this synthetic sig- 
nal are then computed using the first, second, third and 
fourth degree difference statistics {d = 1, 2, 3 and 4) as 
shown in Fig. [2^. We observe that, when using the third- 



degree statistics {d = 3), one has S2 (t) 

.(3) 



C^ 
T^-2 



with 



C2"''' =3.33 that is in good agreement with the theoreti- 
cal value C2 = 2(ci — C2) = 3.3, the classical relationship 
Q = n — 1 being thus well satisfied. When d = 4, we 

obtain Q = 3.35 a value that is consistent with the one 

(2) 
found with d = 3. For d = 2, one obtains Q = 3.23, a 

value that is slightly below the previous one since it be- 
gins to be biased by the smoother part of the signal (two 
times differentiable or more) . Finally, using the usual first- 
degree statistics (d = 1) leads to Q — 2 since the RWC 
signal is differentiable with a steep power spectrum of ex- 
ponent n > 3 (see Table [IJ. 

The estimations of the structure- function exponents Q, 
versus p using the first, second, third and fourth degree 
difference statistics {d = 1, 2, 3 and 4) are presented in 
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Fig. 3: Structure function exponents, (^p, as a function of p 
for tlie synthetic signal (a) and experimental data (b). (o) Q 
computed from the first-degree increments ; (D) Q computed 

from the second-degree differences and fitted by ( ) Q — 

cip-C2pV2 with ci = i.92 and C2 = 0.27; (o) Cp^' and (A) Cp'*' 
computed from the third-degree and fourth-degree differences. 
Solid lines are linear fits of Q : (a) Q = 0.95p, and (b) 
C^'^ = 0.8p. 

Fig. [3^. We observe that Q is a non-linear function of 
p which provides a clear evidence of the intermittent na- 

(3) 

ture of the RWC signal fluctuations. The fact that Q ' 
and Q ' estimates are consistent with the previous ones 



5t(>))/<Tx 



i>-"' (,)/„, 



42) 



>(3) 



.(4) 



{Q ' c^ Q ' ~ Q ') provides further confidence on this di- 
agnosis. It is noteworthy that Q ; Cp and Q estimates 
are in good agreement with the theoretical expectation 
Cp = cip — C2p^/2, which illustrates that the proposed 
framework allows an accurate characterization of the in- 
termittent properties. Finally, as expected for a differen- 
tiable signal, using first-degree increments leads to Q ' is 
a linear function of p with a slope close to 1. This latter 
result thus leads to a misleading conclusion that the RWC 
signal does not present intermittency. This exeinplifies the 
need to adapt the degree of the difference statistics to the 
steepness of the power spectrum. 

Finally, the probability density functions (PDFs) of the 
first-degree increments of the synthetic signal are plotted 
in Fig. [3^ for different time lags r. All the PDFs have 
the same shape independent of the scale r, thus showing 
no intermittency. In contrast, the PDFs of the second- 
degree increments displayed in Fig. |4l3 show a clear evo- 
lution across scales, highlighting the intermittency of the 
signal. This is consistent with the previous structure func- 
tion analysis, and further underlines that high-degree dif- 
ference statistics is needed to test intermittency of steep 
power-spectrum signals. 

Wave turbulence data. We now apply the above pro- 
posed statistical estimator on an experimental signal of 
hydrodynamics surface wave turbulence |20| . A typical 
signal is the temporal evolution of the surface wave ampli- 
tude, ri{t), measured at a given location of the free surface 
of the fluid (see inset of Fig. [T]). Data are recorded from 10 
successive experiments of 300 s each where surface waves 



Fig. 4: PDFs of the first-degree increments 5^ rj for different 
T from 6 to fOO ms (see arrows): (a) Synthetic data, and (c) 
Experimental data. PDFs of the second-degree differences 5^ 77 
for 6 < r < 100 ms: (b) Synthetic data, and (d) Experimental 
data. Dashed-line: Gaussian with zero mean and unit standard 
deviation, a^ axe the rms values. Each curve has been shifted 
vertically for clarity. 

are generated by a wave maker driven by random noise 
forcing in a frequency range 0-6 Hz [50] . Wave heights was 
measured at 1 kHz sampling rate (A = 1 ms) resulting in 
3 X 10^ data points. As for the RWC signal, the initial 
estimation of the power spectrum steepness indicates that 
3 < n < 5. In order to probe for possible intermitent prop- 
erties of this signal with such a steep power spectrum, we 
thus need to use the adapted difference statistics proposed 
in the previous section. 

The power spectrum of 77 (t), estimated using the Fourier 
transform [see Eqs. © and ([7])] of S^ rj, is shown in 
Fig. [1] It displays two frequency ranges with a power- 
law behaviour. In the low-frequency spectrum range 
(~ 7 — 30 Hz) corresponding to the gravity wave turbu- 
lence regime, we observe E{oj) oc a;~ while the high- 
frequency range (~ 30 — 100 Hz) corresponding to the 
capillary regime is characterized by E{uj) oc u!~^'^ j7|. The 
second-order structure functions S2 (t) of rj{t) for d = 1, 
2, 3 and 4 are shown in Fig. [2h as a function of the time 
lag 5 < r < 100 ms. When considering structure func- 
tions 5*2 for d > 2 in the time lag range t < 80 ms 
(corresponding to frequencies above the maximal forcing 
frequency 6 Hz), we roughly observe two different power- 
law scaling behaviours in the gravity regime 15 < r < 65 
ms and the capillary regime 5 < r < 15 ms [a time lag 
T corresponds to a frequency / = 1/(2t)]. We will fo- 
cus below only to the gravity regime since the transition 
between both regimes in Fig. [2b occurs rather smoothly 
which significantly reduce the time lag range available to 
fit the scaling exponent in the capillary regime. As ex- 
plained above, for d = 1, 5*2 is dominated by the signal 
differentiability and does not display both scaling regimes. 
Focusing only on the gravity regime (15 ^ t < 65ms), 
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.(3) 



,(4) 



one consistently finds Q = 3.3 and Q = 3.4 in good 
agreement with the spectral exponent n = 4.3. When us- 
ing d = 2, the smoother transition observed between the 
scaling regimes leads to a slightly underestimated value 

We then look for possible intermittent properties of the 
turbulence wave data in the gravity regime. The evolu- 
tion of Q ' with p is shown in Fig. ^p for d = 1, 2, 3 and 



4. When using the first-degree increments, Q 



(1) 



0.8p 



is a linear function of p. As underlined above, in pres- 
ence of a steep power spectrum, Q, is dominated by the 
differential component of the signal masking possible in- 
terniittency. For d = 2, 3 and 4, we observe a clear non 
linear behaviour of Q, versus p. We note that while Q, 
and Q provide consistent estimates of C,p for all p, Q 
estimates are slightly below these latters as already ob- 
served above for p = 2. Finally, the coherence between 
the estimates of C,p for two successive values of the differ- 

~ Cp ) s-iid with the spectral analysis 



ence degree [Q 

(C: 



(3,4) 



ments are reliable. We can thus conclude these data of 
wave turbulence are intermittent. Fitting of Q with the 
polynomial model C,p ~ cip — C2P^/2 yields ci = 1.9 and 
an intermittency coefficient C2 — 0.27. 

Another way to highlight the wave turbulence intermit- 
tency is to observe a shape deformation of the probability 
density functions (PDFs) of the signal increments with the 
time lag t. The PDFs of the first and second-degree incre- 
ments of the wave amplitudes are respectively plotted in 
Figs.HJ: and|3Ji for 6 < t < 100 ms. When r is increased, 
the PDF's shape of the second-degree increments changes 
continuously up to a nearly Gaussian shape at large r (see 
Fig. [Hi) . This deformation is a direct signature of inter- 
mittency. As predicted above, this intermittency is not 
diagnose when using the first-degree estimator: almost no 
deformation of the PDF shapes of the first-degree incre- 
ments is observed in Fig. [4^. 

In this Letter, we have proposed an easily applicable 
framework based on high-degree difference statistics to 
probe for possible intermittency of a signal with a steep 
power spectrum. We applied it to synthetic data and to 
wave turbulence data. This has led to the observation of 
wave turbulence intermittency |20j . In the same way, it 
can be used on previous existing data notably of magneto- 
hydrodynamic H] or two-dimensional turbulence |B], both 
showing steep power spectra. 
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